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The early Universe is a precious probe of the birth of primordial objects, first star formation 
events and consequent production of photons and heavy elements. Higher-order corrections to 
the cosmological linear perturbation theory predicts the formation of coherent supersonic gaseous 
streaming motions at decoupling time. These bulk flows impact the gas cooling process and de¬ 
termine a cascade effect on the whole baryon evolution. By analytical estimates and N-body 
hydrodynamic al chemistry numerical simulations including atomic and molecular evolution, gas 
cooling, star formation, feedback effects and metal spreading for individual species from differ¬ 
ent stellar populations according to the proper yields and lifetimes, we discuss the role of these 
primordial bulk flows at the end of the dark ages and their detectable impacts during the first 
Gyr in view of the upcoming SKA mission. Early bulk flows can inhibit molecular gas cooling 
capabilities, suppressing star formation, metal spreading and the abundance of small primordial 
galaxies in the infant Universe. This can determine a delay in the re-ionization process and in the 
heating of neutral hydrogen making the observable HI signal during cosmic evolution patchier 
and noisier. The planned SKA mission will represent a major advance over existing instruments, 
since it will be able to probe the effects on HI 21-cm at z ^ 6 — 20 and on molecular line emis¬ 
sions from first collapsing sites at z ^ 20 — 40. Therefore, it will be optimal to address the effects 
of primordial streaming motions on early baryon evolution and to give constraints on structure 
formation in the first Gyr. 
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1. Introduction 

The initial phases of the Universe are extremely interesting for modern Astrophysics, because 
they are supposed to witness the birth of the first stars and galaxies. Features and details of the 
mechanisms leading to this event are still unknown and the properties of primordial generations 
are largely debated (Maio et ah, 2011b; Biffi & Maio, 2013; de Souza et ah, 2013, 2014; Wise 
et ah, 2014). Primordial star formation processes determine the production of photons in different 
energy bands, marking the end of the dark ages and the start of the cosmic dawn. Ongoing feedback 
from stellar evolution impacts the following baryon history and cause the occurrence of the first 
heavy elements in the Universe. Such processes are crucial in the build-up of small primordial 
visible objects, such as gamma-ray-burst host galaxies (Campisi et ah, 2011; Salvaterra et ah, 2013; 
Ghirlanda et ah, 2013; Maio & Barkov, 2014), gaseous damped Lyman-alpha absorption systems 
(Simcoe et ah, 2012; Maio et ah, 2013) and very early faint galaxies (Jeeson-Daniel et ah, 2012; 
Dunlop et ah, 2013; Ouchi et ah, 2013; Maio & Viel, 2014). 

In addition, the emission by atomic (Wouthuysen, 1952; Field, 1958; Bahcall & Ekers, 1969) and 
molecular (Shchekinov, 1986, 1991; Ciardi & Ferrara, 2001; Kamaya & Silk, 2002) lines in the 
pristine gas would be strongly dependent on the thermal state of the high-redshift (z) inter-galactic 
medium. This whole picture is affected by the properties of the initial perturbations that can leave 
their imprint in the infancy of primordial structures. Indeed, supersonic streaming motions are 
expected to be originated as a result of the relative velocities between dark matter and baryons 
following recombination, at z ~ 10^ (Tseliakhovich & Hirata, 2010). They are due to non-linear 
corrections to cosmological perturbation theory that predicts coherent gas motions on Mpc scales 
with typical velocities with a rms value of ~ 30km/s. These bulk flows could be able to influence 
cosmic gas emission (McQuinn & O’Leary, 2012; Pritchard & Loeb, 2012) and star formation 
events at z ~ 10 — 30 (Maio et ah, 2011a; Stacy et ah, 2011; Greif et ah, 2011; Maio & lannuzzi, 
2011; Fialkov, 2014), with possible consequences on the kinetic SZ effect and on the determination 
of the temperature-density equation of state (Dalai et ah, 2010). Although direct observations 
are still very challenging and outside the possibilities of current facilities, they could be probed 
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by the future Square Kilometre Array (SKA) (Carilli & Rawlings, 2004; Dewdney et al., 2013) 
that will deteet radiative emissions by cosmic gas at redshift z ~ 6 — 30. At those epochs the 
Universe is largely pristine, dominated by H atoms that can form H 2 and HD and lead to gas cooling 
and collapse. Molecular-rich gas in dense regions is likely to emit through collisional excitations 
around ~ 533 GHz and such features can be used to identify the very first star forming sites. The 
following star formation processes would locally inject large amounts of entropy (e.g. Maio et al., 
2013) and would produce energetic photons that will boost HI 21-cm (1.4 GHz) emission via Ly- 
a pumping, allowing the gas to be seen against the cosmic microwave background (Wouthuysen, 
1952; Field, 1958; Madau et al., 1997; Wyithe & Loeb, 2003; Ciardi & Ferrara, 2001). The resulting 
temperature contrast in the optically thin limit depends on the optical depth, T, which is a function 
of the neutral H fraction. Therefore, 21-cm analysis is a very suitable probe of the ages before re¬ 
ionization and can help place constraints on cosmic star formation, early galaxy evolution and on 
the implications of primordial bulk flows (Dalai et al., 2010; McQuinn & O’Leary, 2012; Visbal et 
al., 2012; Mellema et al., 2013). Given the fact that SKA will observe in different frequency ranges 
- between 50 MHz (SKAl-LOW) and 14 GHz (SKAl-MID) in Phase 1 and up to 24 GHz at 10 
times larger sensitivity in Phase 2 - these radiative emissions will result within its instrumentational 
potential and, if detected, will represent a strong advancement for our theoretical understanding of 
primordial structures. 

In the following, the impacts of baryon streaming motions for early gas emissions and structure 
formation will be outlined and the possibilities for the different frequency ranges of SKA will be 
discussed in more detail. 

2. Streaming motions at the of the dark ages and perspectives for SKA 

Streaming motions affect in the first place the primordial gas evolution, delaying its collapse 
and the subsequent star formation history. This is shown by ad hoc numerical simulations (Maio et 
al., 2011a) adopting a ACDM model with present-day geometrical parameters for total-matter. A, 
dark-matter and baryon density given by: = 0.3, flo,A = 0.7, Do,dm = 0.26, Do,b = 0.04, and 

an expansion parameter at the present Hq = 100/ikm/s/Mpc with h = 0.7. A simple way to take 
into account primordial bulk velocities in numerical simulations is to assume a velocity shift along 
the e.g. x-direction, Vbx, given uniformly to each gas particle at initial time, according to magni¬ 
tudes of Vbx = 0, 30, 60 km/s at z ~ 1020 - see (Naoz et al., 2011) for further studies. In this way 
we explore various streaming velocity scenarios by means of the numerical code Gadget (Springel, 
2005), with subsequent modifications including: non-equilibrium cosmic chemistry (for e^, H, 
H+, H , He, He+, He++, H 2 , H^, D, D+, HD, and HeH+); stellar evolution of stars for different 
masses, metal yields (He, C, O, Mg, S, Si, Fe, N, etc.) and lifetimes (Maio et al., 2007; Tornatore 
et al., 2007; Maio et al., 2010); star formation led by H, He and H 2 radiative losses in pristine envi¬ 
ronments and by resonant or fine-structure metal transitions (such as CII 157.7 /rm. Sill 34.8 /rm, 
01 145.5 and 63.18 /rm. Fell 87.41, 51.28, 35.35, 25.99 /rm) in enriched environments (Maio et 
al., 2007). Subsequent supernova and wind feedback are responsible for spreading heavy elements 
in the surrounding gas and for determining the transition from the primordial pristine (population 
III) star formation regime to the metal-enriched solar-like (population II-I) regime when a critical 
metallicity of lO^"* Z© is reached - see further details in Maio et al. (2010). Cosmic structures at 
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Figure 1: Expected differential (left column) and cumulative (right column) distributions for gas clouds at 
redshift z = 23 (top) and 19 (bottom). Data refer to simulations of ^1 Mpc-size box run with different (7g 
and Vbx'- Cg = 0.9 and v^x = Okm/s (solid line and crosses); CTg = 0.9 and Vhx = 60km/s (dotted line and 
triangles); CJg = 0.8 and V],x = Okm/s (dashed line and asterisks); CJg = 0.8 and Vhx = 60km/s (dot-dashed 
line and squares). 


different redshifts are identified by means of a friends-of-friends algorithm with a linking length of 
20% the mean interparticle separation. For each object gas, dark and stellar properties are stored. 

2.1 Implications for gas evolution 

Fig. H] displays the differential (left) and cumulative (right) distribution of gas clouds at redshift 
z = 23 (top) and 19 (bottom), corresponding to cosmic times of about 0.1 and 0.2 Gyr, respectively, 
in different parameter scenarios. They refer to four simulations run with initial conditions generated 
for a 0.5 Mpc//i-side box, sampling dark-matter and gaseous-matter fields wifh 2 x 320^ parficles 
and corresponding fo a gas resolution of abouf 40 Mq. A primordial specfral slope of n = 1 is 
assumed and fhe normalizafion is fixed by imposing ag = 0.9. In order fo fesf degeneracies wifh 
specfral paramefers, resulfs for a case of ag = 0.8 are overploffed, as well. In general, sfreaming 
motions confribufe in shaping fhe gaseous confenf and in lowering fhe typical clump masses and 
abundances. This is mosfly evidenf af higher redshifl when dark-mafler haloes are smaller and fhe 
addifional kinefic energy coming from fhe bulk flows make more difficull fo frap fhe gas in fhe 
dark-mafler potential wells. Panels referring fo z = 23 highlighf fhis effecl very clearly. The differ- 
enfial disfribufion shows fhal, for any choice of ag, the v^x = Okm/s case predicts larger and more 
abundant gaseous structures, while the v^x = 60km/s scenario predicts typically lower-mass struc¬ 
tures. This is due to the effects of bulk flows fhal cause gas evacualion from primordial objecls. If 
is worlh nofing fhal fhe effecls of sfreaming velocifies is comparable fo a significanl decrease of ag 
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from 0.9 to 0.8. The overall changes of cosmic abundances are better displayed by the cumulative 
distribution at the same redshift. In the Vbx = Okm/s case variations in the high-mass tail of the 
number densities of up to 1 dex with respect to the limiting case of Vbx = 60km/s are reported. In 
the lower-mass end suppression of a factor of ~ 2 is visible. The slight excess of small clumps in 
the Vbx > Okm/s models is due to residual gas that collapses later with respect to the Vbx = Okm/s 
case. At later times, it is still possible to see similar trends, as shown in the z = 19 panels. There 
differences between different-Vb^r models decrease slightly due to the simultaneous increase of typ¬ 
ical dark-matter masses and the redshift decay of bulk flows. Thus, gas in-fall and collapse are less 
hindered. This leads to typical differences of a factor < 2 in the gas cloud differential distributions 
and of a factor of a few in the gas cloud cumulative distributions. The statistical behaviour shows a 
converging evolution of the gas cloud masses and abundances at later times (z ~ 10 — 20). At these 
epochs gas collapse is only marginally affected by primordial streaming motions (at ~ 10% level) 
and a major role is played by cosmological star formation and environments (see next). These 
trends suggest that the main effect on gas cloud evolution is a delay of the collapse redshift of Az ~ 
a few (corresponding to tens of Myr at the epoch of interest). 

Since in these epochs there are no significant star formation events, yet, and forming structures 
are characterized by abundant pristine molecular gas, SKA could be able to detect H 2 and HD 
signatures from first collapsing sites. E.g., H 2 7 = 2 — 0 rotational transition and HD 7 = 4 — 3 
rotational transition have rest-frame frequencies of 533.179 GHz and 533.388 GHz, respectively. 
During Phase 1 SKA 1-MID will have frequency bands up to ~ 14 GHz, hence observable (red- 
shifted) signal from rare dense collapsing gas at very high redshift (around z ~ 35 — 40) will fall 
in the covered regime. However, the SKA operational range will reach 24 GHz in phase 2, and 
then SKA frequency bands will cover the frequencies of molecular lines emitted by cold gas clouds 
down to z > 20. The emissivity of such lines should be detactable. Indeed, H 2 7 = 2 — 0 emission 
is expected to be at level of ~ 2 x 10'^ Jy per each collapsing clump, while HD 7 = 4 — 3 line is 
expected to have a lower rate of ~ 10^^ Jy (Kamaya & Silk, 2003). Roughly speaking, our theoret¬ 
ical estimates in Fig. |] suggest a number of collapsing clumps ranging from some tens up to some 
thousands per Mpc^ at z ~ 23 — 19, hence the total emissivity by a ~ 1 Mpc^ patch (accounting 
for an area of r? ~ 0.33deg^ at z ~ 19)^ should range from ~ 10 up to ~ 1000 times 2 x 10~’Jy, 
i.e. ~ 2 — 200/iJy. Assuming field of view and sensitivity for SKAl-MID (A = 0.49deg^ and 
63/xJy — hr~'/^, respectively) the detectable emission turns out to be A/rJ = 0.49/0.33 ~ 1.5 times 
higher, i.e. around ~ 3 — 300/^Jy. Therefore, the effects could be assessed by a-few-hours SKAl 
observations of gas molecular lines. The signal should be more easily detected by SKA2-LOW that 
is going to have a 10-times better sensitivity in the 350 — 24000MHz frequency band. 

2.2 Implications for star formation 

After the condensation of first gas clumps, star formation processes start taking place and 
radiation from stars act on the thermodynamical state of the surrounding medium. 

The top panel of Fig. ^ displays the effects of streaming motions on the whole cosmological star 
formation rate density and on the population III (pop III) contribution. In all the cases there is 

*The (comoving) distance to z ~ 19 is about lOGpc, the corresponding angular separation of 1 (comoving) Mpc re¬ 
sults to be 10^^rad ~ 0.573deg and hence 1 Mpc^ patch covers 0.33deg^. This is comparable to SKAl-MID resolution. 
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Figure 2: Top panel: Redshift evolution of the total cosmological star formation rate density and of the 
population III (pop III) star formation rate density for Vb.j = Okm/s (solid lines), Vbx = 30km/s (dotted lines), 
Vhx = 60km/s (dashed lines). Bottom panel: Contribution of the pop III star formation regime (SFRni) to 
the cosmological star formation rate (SFRtot) for the different Vb^r cases. 


a cosmological effect on the trend and magnitude of the star formation activity. This suffers the 
delays in the gas collapse caused by streaming motions and shows a retarded onset for different Vbx 
values, up to Az ~ 3. For the Vb^ = Okm/s case the onset takes place at z ~ 20, when the age of 
the Universe is about 0.17 Gyr, while for the cases including a velocity shift star formation kicks 
in at z — 19 and z — 17.6 for Vbx = 30km/s and Vb,: = 60km/s, respectively. The corresponding 
cosmological times are 0.19 Gyr, and 0.21 Gyr accounting for time delays of some ten Myr with 
respect to the v\,x = Okm/s case. Values for the star formation rate density differ at very early 
times of more than one order of magnitude and global trends show a converging behaviour only 
below z ~ 14, catching up at z < 10. The reason for this is related to gas cloud evolution that 
results delayed because of the higher kinetic energy gained at recombination. In larger-Vbjc models, 
primordial mini-haloes can retain less material (see Fig. |]), the hosted gas condenses more slowly, 
molecule formation gets hindered, the consequent cooling is less efficient and the resulting star 
formation activity is delayed. This impacts mainly primordial haloes with masses 10"^ - 10^ Mq, 
as their dimensions are comparable to or smaller than the baryon Jeans length (Tseliakhovich & 
Hirata, 2010). Therefore, gas cannot fragment within the halo and partially or entirely flows out 
of it, as larger dark-matter potential wells would be needed to completely retain it. Basically, 
these processes start from small scales and cascade over the growth of larger objects suppressing 
gas condensation and subsequent star formation already before reionization. Photon production, 
feedback effects and the consequent occurrence of the first heavy elements in the Universe will be 
delayed, as well. This is better seen in the lower panel of Fig. ^ that highlights the contribution 
from the pristine pop III regime (whose evolution is dictated by ongoing metal production) to the 
cosmic star formation rate. It is evident that for larger Vbx the transition from the primordial pop 
III era to the metal-enriched pop II-I one takes place at later times due to the impacts of early bulk 
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flows. At later times (z < 10) the trends eonverge and the imprints of the initial dynamieal patterns 
are more diffieult to reeognize. 

2.3 Implications for 21-cm signal 

A valuable observable to investigate the eosmologieal response of baryons at the end of the 
dark ages (z ~ 6 — 20) is 21-em signal from HI. This falls in the frequeney range eovered by SKAl- 
LOW (50 — 350 MHz), that eould possibly map at a 3—a level neutral gas from the dark ages and 
eosmie dawn until full re-ionization at z ~ 6 (Mellema et ah, 2013). 

The expeeted signal is usually deseribed in terms of differential brightness temperature, 8T},, the 
extent of whieh depends mainly on the ionization and thermodynamieal state of the (partially) 
neutral gas. The effeet on star formation eaused by primordial streaming motions will impaet 
8T}, in different ways. In partieular, by delaying the onset of star formation streaming motions 
would also shift the beginning of the re-ionization proeess towards lower redshift (i.e. higher 
frequeneies), while by suppressing or redueing star formation within small mass halos they would 
ehange the topology of re-ionization. As a eonsequenee, an impaet on the 21-em signal is also 
expeeted. Semi-analytieal arguments extrapolated to ~ 100Mpe seales (Dalai et ah, 2010; Visbal 
et ah, 2012; MeQuinn & O’Leary, 2012; Ali-Haimoud et ah, 2014) suggest Vb^ implieations on the 
21-em signal from eomoving distanees of the order of ~ 100kpe (Visbal et ah, 2012; MeQuinn & 
O’Leary, 2012) down to ~ lOkpe (MeQuinn & O’Leary, 2012; Ali-Haimoud et ah, 2014), with an 
overall effeet on 8T\^ ranging from ~ 2mK (Dalai et ah, 2010) to ~ lOmK (Visbal et ah, 2012). On 
the other hand, implieations from strueture formation shoeks ean eontribute at < 10% level, while 
enhaneements to the ionization fraetion due to X-rays, albeit very uneertain, seem to be negligible 
or almost irrelevant (e.g. MeQuinn & O’Leary, 2012). 

We finally note that the seenario eommonly adopted for these estimates relies on standard eold dark 
matter, however, eomparable or even larger suppressions of eosmie strueture evolution might arise 
in presenee of alternative baekground models, sueh as warm dark matter (see e.g. Maio & Viel, 
2014, and referenees therein), non-Gaussian matter distributions (Maio & lannuzzi, 2011; Maio et 
ah, 2012; Maio & Khoehfar, 2012; Paee & Maio, 2014) or early dark energy (Maio et ah, 2006). 

3. Conclusions 

Higher-order eorreetions (Tseliakhovieh & Hirata, 2010) to the eosmologieal linear perturba¬ 
tion theory prediet the formation of eoherent supersonie gaseous streaming motions at deeoupling. 
These are relevant for primordial gas eondensation, the birth of the first stars and the build-up of the 
lowest-mass objeets and proto-galaxies. In view of the upeoming SKA, we have diseussed the im¬ 
plieations of sueh bulk flows via numerieal N-body hydrodynamieal ehemistry simulations (Maio 
et ah, 2011a) of early strueture formation taking into aeeount non-equilibrium eosmie ehemistry, 
gas eooling (Maio et ah, 2007), stellar evolution of stars for different masses, metal yields (He, C, 
O, Mg, S, Si, Fe, N, ete.) and lifetimes (Maio et ah, 2007; Tornatore et ah, 2007) and the transition 
(Maio et ah, 2010) from the primordial pristine pop III star formation regime to the metal-enriehed 
solar-like pop ILL Primordial gas streaming motions ean suppress gas eollapse, with eonsequent 
delay of eosmologieal star formation, feedbaek effeets, metal spreading and eosmie re-ionization. 
Their impaets on gas elouds are partieularly evident at z ~ 15 — 30, when number densities and 
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typical gas content in primordial haloes result lowered of up to a factor of 10. Hence, they af¬ 
fect molecular H 2 and/or HD emissions from overdense gaseous regions that will be observable by 
SKAl-MID. Subsequently, bulk flows delay star formation, metal pollution and feedback mech¬ 
anisms of Az ~ a few at z ~ 10 — 20, influencing the cosmic HI thermal and chemical state, and 
making the 21-cm signal at the end of the dark ages patchier and noisier. The next SKA mission 
will represent a major advance over existing instruments, since it will be able to detect HI 21-cm 
(e.g. SKAl-LOW at 50 — 350MHz) and molecular-line (e.g. SKAl-MID up to 14 —24GHz in 
phase 1-2) emissions from early star forming sites at z ~ 6 — 40 (SKA2-LOW). Therefore, it will 
be possible to address the effects of primordial streaming motions on early baryon evolution and to 
give constraints on structure formation in the first Gyr. 
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